Monte Carlo Simulations of Short-time Critical 
Dynamics with a Conserved Quantity 

B. Zheng^'^ and H.J. Luo^ 

^FB Physik, Universitdt Halle, D-06099 Halle, Germany 
'^Institute of Theoretical Physics, Academia Sinica, 100080 Beijing, P.R. China 
^FB Physik, Universitdt Siegen, D~57068 Siegen, Germany 

Abstract 

With Monte Carlo simulations, we investigate short-time critical dynamics 
of the three-dimensional anti-ferromagnetic Ising model with a globally con- 
served magnetization rUs (not the order parameter). From the power law 
behavior of the staggered magnetization (the order parameter), its second 
moment and the auto-correlation, we determine all static and dynamic crit- 
ical exponents as well as the critical temperature. The universality class of 
rus = is the same as that without a conserved quantity, but the universality 
class of non-zero rUg is different. 
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I. INTRODUCTION 



About ten years ago, in a pioneering work by Janssen, Schaub and Schmittmann, short- 
time universal scaling behavior of non-equilibrium critical dynamics starting from a high 
temperature state was systematically explored with renormalization group methods In 
a few years, extension of the calculations to different critical dynamics was carried out by 
Oerding and Janssen in a series of works 0-^. Some evidences for the short-time dynamic 
scaling were also observed in Monte Carlo simulations Meanwhile, it was found that 

the power law decay of the magnetization starting from a completely ordered state emerges 
already in relatively early times, e.g., see [^, and it can be used to estimate the dynamic 
exponent z PJTO[|. In recent years short-time critical dynamics has been systematically in- 
vestigated with Monte Carlo simulations JTT|-[T3[]. Simulations have extended from regular 



classical spin models |T^J16| -[T9| to statistical systems with quenched disorder p0|-[2^, quan- 
tum spin systems lattice gauge theories the hard-disk model P3|PB| and dynamic 
systems without detailed balance p7|-|3T[] . The references given here are only a part of recent 
ones and not complete. A relatively complete list of the relevant references before 1998 can 
be found in a recent review, Ref. All numerical and analytical results confirm that 

there exists a rather general dynamic scaling form in critical dynamic systems already in 
the short-time regime of dynamic evolution. 

Traditionally, it is believed that the short-time behavior of dynamic systems depends 
essentially on microscopic details. To understand short-time universal behavior, one should 
distinguish two different time scales, the macroscopic and microscopic time scales. The 
short-time dynamic scaling form emerges only after a time scale tmic, which is sufficiently 
large in microscopic sense, but still very small in macroscopic sense, tmic is a time that the 
system needs to sweep away the effect of microscopic details. In Monte Carlo simulations, 
for example, if a sweep over all spins on a lattice is considered to be a microscopic time 



unit, tmic is usually the order of 10 to 100 Monte Carlo time steps [|1J]. For the simple Ising 



and Potts models with only nearest neighbor interactions, sometimes tmic is negligible small 
(e.g., one Monte Carlo time step or two). Compared with typical macroscopic time scales 
characterized by t~^^ or around the critical point, tmic is indeed very small. 

The physical origin of the short-time dynamic scaling is the divergent correlation time 
near a critical point. The divergent correlation time induces a memory effect, and the 
memory effect can be described by a scaling form. How does this scaling form look like? 
As pointed out by Janssen, Schaub and Schmittmann |]l|, one should introduce new critical 
exponents to describe the dependence of the scaling behavior on the macroscopic initial 
conditions. Actually, for arbitrary initial conditions even a characteristic function is needed 
|5^,|55|]. Furthermore, to describe the scaling behavior of some special dynamic observables. 



we also have to introduce new exponents, for example, the persistence exponent [pHp5 



The short-time dynamic scaling form is not only conceptually interesting but also prac- 
tically important. It provides new techniques for the measurements of both dynamic and 
static critical exponents as well as the critical temperature, for a review, see Ref. |r5]. Since 
the measurements are carried out in the short-time regime of dynamic evolution, the dy- 
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namic approach does not suffer from critical slowing down. [| Averaging is over both initial 
configurations and the random forces. This is very different from the time averaging in the 
measurements in equilibrium. Compared with those non-local methods which are invented 
to overcome the critical slowing down in Monte Carlo simulations in equilibrium, e.g., the 
cluster algorithms, the dynamic approach does study the original local dynamics and can 
be applied to disordered systems. 

In fact, the dynamic scaling in non-equilibrium critical dynamics is not an exclusive 
phenomenon in the world. In many other non-equilibrium dynamic processes, universal 
or quasi-universal scaling behavior has been also observed. An example is phase ordering 
dynamics In this real equilibrium state and a macroscopic time scale like 

T~'^^ do not exist. Therefore, 'short-time' is not addressed. Another example is aging in 
complex systems such as glasses and spin glasses. Aging is just a kind of scaling or quasi- 
scaling. Dynamic scaling behavior around a spin glass transition P,|37 H39|j2l[| is very similar 
to that around a standard critical point. For example, the experimental measurements of 
the remanent magnetization in spin glasses support not only the power law scaling behavior 
but also the scaling relations between the exponents 

After such a slightly lengthy review, we come to what we are going to do in this paper. 
Even though renormalization group calculations have been extended to different critical 
dynamics, up to now, Monte Carlo simulations are limited in dynamics of model A El 



Actually, due to severe critical slowing down, Monte Carlo simulations for long-time behav- 
ior of critical dynamics beyond model A are also in the preliminary stage. The short-time 
dynamic approach provides powerful methods for numerical measurements of not only dy- 
namic exponents, but also static exponents as well as the critical temperature. Simulations 
of dynamic systems beyond model A in equilibrium are usually difficult. Therefore, it is 
important to investigate short-time scaling behavior of these dynamic systems with Monte 
Carlo methods. 

Dynamics of model A is a kind of relaxational dynamics without relevant conserved 
quantities If we consider only dynamic relaxational processes, dynamics with different 
relevant conserved quantities is classified into model B, C and D. In this paper, taking the 
three-dimensional anti-ferromagnetic Ising model for example, with Monte Carlo methods 
we study dynamics with a conserved quantity, which is not an order parameter but coupling 
to the order parameter. According to Ref. this is called model C. In the next section, 
we introduce the model and analyze the power law scaling behavior in the short-time regime. 
In section 3, we present the numerical results. Finally come the concluding remarks. 



II. SHORT-TIME DYNAMIC SCALING 



In the literature, sometimes it is stated that the short-time dynamic approach can eliminate 
critical slowing down. Rigorously speaking, this statement is not correct. Critical slowing down 
always exists, at least for local dynamics, and it is the physical origin of the dynamic scaling. 
But in the short-time dynamic approach we do not have the problem of generating independent 
configurations. Therefore, we do not suffer from critical slowing down. 
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A. The model 



We consider an antiferromagnetic Ising model on a three-dimensional cubic lattice. The 
Hamiltonian of the model is 

- H/kT = -Kj2 S^S„ (1) 

<ij> 

where Si is an Ising spin and the sum is over nearest neighbors. For dynamics of model 
A and in equilibrium, on a cubic lattice the antiferromagnetic Ising model is equivalent to 
the (ferromagnetic) Ising model. Now, we keep the magnetization as a constant in dynamic 
evolution. Here the order parameter is not the magnetization but the staggered magneti- 
zation. We will denote the magnetization by and the staggered magnetization by M{t). 
For equilibrium states, it is expected that for = 0, the universality class is the same as 
that of the standard Ising model without a conserved magnetization. For a non-zero m, 



according to Ref. the critical exponent z/ is different from that of = 0.0, but it does 
not depend on the value of the non-zero m^. For critical dynamics, the dynamic exponent 
z and the exponent 6 for = are the same as those of the standard Ising model, but for 
non-zero they are different. Again, z and 6 do not depend on the value of the non-zero 
mg. These conclusions are drawn from renormalization group calculations based on the 0^ 



theory coupling to a conserved current ||4ll,p[|. It is interesting to compared these results 
with Monte Carlo simulations. 

In this paper, we study only dynamic relaxation starting from disordered states with a 
zero or small initial staggered magnetization ttiq. We measure how the staggered magneti- 
zation and its second moment as well as the auto-correlation evolve in the dynamic process. 
The heat-bath algorithm is used in simulations. In order to keep constant, in a flip 

we simply exchange the values of two spins. Here we should point out that if our model is 
not the anti-ferromagnetic but the ferromagnetic Ising model, the (conserved) magnetization 



rus is the order parameter and then the dynamics for = belongs to model B |^ . 

In dynamics without any relevant conserved quantities, i.e., dynamics of model A, up- 
dating schemes are irrelevant in the sense of universality. Different updating schemes lead to 
the same critical exponents, either in an equilibrium or a non-equilibrium state. However, 
for dynamics with relevant conserved quantities, it is different, at least for non-equilibrium 
short-time behavior. In a recent work by one of the authors [Q, for example, the Monte 
Carlo dynamics with a conserved order parameter for the two-dimensional Ising model is 
investigated. It is found that if in a flip we exchange the values of two neighboring spins, the 
dynamics is very slow. However, if we release the condition such that the order parameter 
is not locally conserved but only globally conserved (for example, in a flip we exchange the 
values of two randomly separated spins), the dynamic exponent z is much smaller. For 
the dynamics studied in this paper, i.e., with a conserved quantity which is not the order 
parameter, the situation is somewhat similar. If in a flip we exchange the values of two 
neighboring spins, dynamic evolution is somewhat slow. Therefore, as a first approach, we 
decide to keep the magnetization only globally conserved. In a flip we exchange the val- 
ues of two randomly separated spins. Except for the updating schemes, all computational 
techniques adopted here are the same as those in the simulations of model A WM. 
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B. Short-time dynamic scaling 



For critical dynamic systems, traditionally it is believed that universal scaling behavior 
exists only in the long-time regime of dynamic evolution. However, in recent years it is 
discovered that starting from macroscopic initial states, universal scaling behavior emerges 
already in the macroscopic short-time regime of dynamic processes after a microscopic time 
scale tmic |[l|,|5|J^,p,12,14]. A typical example is that a magnetic system initially in a high 



temperature state with a small initial order parameter mg, is suddenly quenched to the 
critical temperature or nearby (without external magnetic field) and then released to 
dynamic evolution of model A [^. A generalized dynamic scaling form can be written 
down, for example, for the kih moment of the order parameter, 

M^^^ {t,T,L,mo) = fo-'^/^/^^Cfc) ^^~z^^ ^i/v^^ h-^L,}f°mo). (2) 

Here t is the time variable, r is the reduced temperature, L is the lattice size. (3 and u 
are standard static exponents and z is the dynamic exponent. Important is that a new 
independent exponent Xq is introduced to describe the scaling behavior of the initial order 
parameter itlq. If the scaling form above is valid, all relevant exponents can be extracted 
from the short-time behavior of relevant observables. 

For the dynamic system with a conserved quantity discussed in the last subsection, we 
assume that a scaling form like Eq. holds also. In this case, r = {K — Kc)/Kc and 
Kf. depends on the conserved magnetization m^, i.e., there is a critical line = Kc{ms). 
In principle, there are two possibilities for the critical exponents. In the case of strong 
universality, the value of an exponent for non-zero can be different from that for zero 
rris, but does not depend on the value of the non-zero m^. In the case of weak universality, 
an exponent may vary its value continuously along a critical line. For the (j)^ theory coupling 
to a conserved current in Ref. it is the case of strong universality. 



Neglecting the finite size effect and noting that mo is small, we expand the right hand 
side of Eq. and take only the first non-zero (linear) term under the condition of a small 
t. At the initial stage of time evolution, the staggered magnetization around the critical 
temperature behaves as, 

M{t, r, mo) ~ mo F{t^''''T), 6 = {xq - f3/iy)/z. (3) 

At the exact critical point, r = 0, M{t) obeys a power law ~ t^. Numerical results and 
analytical calculations have revealed that the exponent 6 is positive for almost all systems, 
i.e., the order parameter undergoes an initial increase. This makes the short-time behavior 
very prominent. The physical mechanism for this increase has not been very clear. At 
least the mean-field effect or symmetry breaking is not very relevant. Slightly away from the 
critical point, the power law behavior is modified by the scaling function F^t^^'^^r). This fact 
allows us to locate the critical temperature and meanwhile to measure the critical exponent 

e. 

Differentiation of Eq. (^ leads to 

drM{t, r, mo) 1^=0 = ^o 9,,F(r') |.'=o, Q = l/uz + 9. (4) 
From this power law behavior, one can determine the critical exponent 1/vz. 
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Taking into account that the non-equihbrium spatial correlation length (~ t^/^) is small 
at the initial stage of time evolution, from Eq. (|^) it can be derived that the second moment 
at Tc subjects to a finite size scaling 



M^^\t, L) L~'^ ty , y = {d-2l3/v)/z. 



(5) 



For simplicity, here we have set mo = 0. Another interesting dynamic observable is the 
auto-correlation 



A{t) 



(6) 



At the critical temperature Tc, A{t) decays by a power law [^] 



A{t) ~ r 



A 



-9. 



(7) 



It is interesting that even though we have set mo = 0, 6^ (i.e., xq) still enters the auto- 
correlation. This is because xq is actually the scaling dimension of the local order parameter. 

From Eqs. (|^) and (0), we are able to estimate the static exponent l3/v and the dynamic 
exponent z. Then we complete the measurements of all the exponents and the critical 
temperature. For details of the above scaling analysis and a more systematic extension, 
readers are referred to Refs. p|,^ and the recent review article Ref. WA . 



III. MONTE CARLO SIMULATIONS 
A. = 

For the anti-ferromagnetic Ising model, in equilibrium the averaged magnetization is 
zero. If we notice that the order parameter is the staggered magnetization, the thermody- 
namic fluctuation of the magnetization goes to zero in the thermodynamic limit (infinite 
volume). In non-equilibrium dynamic processes without any conserved quantities, if one 
starts from initial states with a zero magnetization m^, the magnetization will remain zero. 
This is the case for usual macroscopic initial states, e.g., the high temperature and low tem- 
perature states. The fluctuation of the magnetization in the dynamic processes is also zero 
in thermodynamic limit. Therefore, one expects that the equilibrium state of the dynamic 
process with a conserved magnetization m^ = is in a same universality class of the model 
without a conserved quantity, and so is the dynamic universality class as well. The critical 
temperature is also not changed. 

We perform simulations for m^ = in order to confirm the above expectation and to 
demonstrate our computational techniques. The heat-bath algorithm is adopted in simu- 
lations since it is faster than the Metropolis algorithm in the short-time regime. Lattice 
sizes are L = 64 and 128. Within our maximum updating times, no finite size effects are 
observed. We remind the readers here that it is an advantage of the short-time dynamic 
approach that the finite size effects are easily controlled due to the small non-equilibrium 
spatial correlation length (~ t^^^). 
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The critical temperature has been measured rather accurately in simulations in equi- 
librium (for the system without a conserved quantity). However, we also present data to 
determine it from short-time dynamics. With the short-time dynamic approach, in princi- 
ple, any observables which are sensitive to the temperature around the critical regime can 
be used for the determination of the critical temperature. 

In Fig. ID, time evolution of the staggered magnetization M{t) is displayed with solid lines 
in log-log scale. The lattice size is L = 128, and 2500 samples of initial configurations have 
been used for averaging. In order to see the possible 'finite rriQ effects', we have performed 
the simulations with mo = 0.01 and 0.02. The coupling constant K is 0.22065, 0.22165 and 
0.22265 (from below) for mo = 0.01, and 0.22115, 0.22165 and 0.22215 for mo = 0.02. In 
principle, the closer i^'s are used, the more accurate Kc and the critical exponent l/uz can 
be obtained. However, we suffer from large statistical fluctuation if the K^s are too close to 
each other. We use different i^'s for mo = 0.01 and mo = 0.02 just to study the possible 
systematic errors. 

To make use of Eq. to locate the critical point, we first interpolate M(t) quadratically 
to any K around the three simulated ones, then search for a K which gives the best power 
law behavior for the curve. This K is Kc. The two dashed lines in Fig. |l] are the curves 
with the best power law behavior for mo = 0.01 and mo = 0.02. The corresponding K^s 
are = 0.22169(9) and 0.22163(5) for mo = 0.01 and mo = 0.02 respectively. These two 
values are consistent each other and in agreement with that of the Ising model measured 
in equilibrium (without any conserved quantity), which is reported to be around 0.22165 
P^ , |45[| . We perform our measurements in a time interval [tmic,400]. The results are rather 
stable when we take tmic bigger than 10. 

With Kc in hand, we measure the critical exponent 6 = 0.108(5) and 0.100(5) for mo = 
0.01 and mo = 0.02. They are consistent with 6 = 0.104(3) in Ref. [|TB[ and = 0.108(2) 
in Ref. |^ for the Ising model without a conserved quantity. Rigorously speaking, the 
exponent 9 is defined in the limit of mo = 0. We should extrapolate 9 to mo = 0. However, 
since here there are some statistical errors, and especially the errors induced by the errors 
of the critical point KcS, it is not meaningful to do so. To reduce the errors induced by Kc, 
we could take Kc = 0.22165 as input. Then we obtain 9 = 0.106(6) by extrapolating the 
results to mo = 0. In this paper, we simply consider 9 = 0.108(5) measured from mo = 0.01 
as our final value of 9. 

In Fig. 0, time evolution of the second moment of the staggered magnetization is dis- 
played with solid lines in log-log scale. Here the initial staggered magnetization mo has been 
set to zero. From below, K = 0.22065, 0.22165 and 0.22265. The lattice size is L = 64 and 
the number of samples for averaging is 3200. Extra simulations for L = 128 show that the 
finite size effect for L = 64 is already negligible small. From the figure, we see that the sec- 
ond moment is apparently less sensitive to K, compared with the staggered magnetization 
itself. However, we still can locate the critical Kc from the data. Measuring in a time interval 
[tmic = 20, 700], Kc = 0.22164(22). The curve corresponding to Kc is shown by a dashed line 
in the figure. The error of Kc here is bigger than that from the staggered magnetization. To 
reduce the error, simply increasing the samples for averaging is not enough, one must have 
longer maximum updating times and therefore larger lattices. These need much CPU times 
since the second moment is not self-averaging. To estimate the exponent y = {d — 2(3/i>)/z, 
for simplicity, we take Kc = 0.22165 as input. The result is y = 0.965(11). 
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In the simulations for Fig. ^, we have also measured the auto-correlation function A{t). 
The results are shown in Fig. |. From below, K = 0.22065, 0.22165 and 0.22265. It is 
impossible to locate the critical point from the auto-correlation. We measure the exponent 
A = 1.36(3) at = 0.22165 in a time interval [20,100]. When t is bigger than 100, the 
fluctuation becomes large. 

In order to determine the critical exponent z/, we need drM{t, r)|^=o- From the data for 
Fig. |I|, we can calculate this derivative approximately. The results are shown in Fig. ^ with 
solid lines in log-log scale. From below, the initial staggered magnetization is mg = 0.01 
and 0.02. It is known that the power law behavior of the derivative of M{t) in shorter 
times is less clean than that of M{t) itself |[T^. Therefore, we measure the slopes in a time 
interval [70,400] and obtain Cd = l/uz + 6 = 0.867(10) and 0.901(22) for mo = 0.01 and 
0.02 respectively. The value from mo = 0.01 is more reliable. 

In Table |, we list all the exponents we have measured. Taking the exponent 6 as input, 
from A we estimate the dynamic exponent z. Then from y and q we calculate the exponents 
2(3 /u and u. For comparison, available results for dynamics without a conserved quantity 
P6|| , i.e., the so-called model A dynamics, are also given, where the critical point Kc and the 
exponent u are measured from another dynamic process starting from a completely ordered 
state. Within statistical errors, all the exponents for dynamics with and without a conserved 
quantity agree well each other. The static exponents and are also consistent with those 
of the standard Ising model measured from simulations in equilibrium with the cluster 



algorithms [^]. Here we should make a comment. Even for dynamics of model A, the static 
exponents and Kc measured from the short-time dynamics are still not as accurate as those 
obtained with the cluster algorithms. One reason is that the best effort for the short-time 
dynamic approach has not been made. Another reason is that in the short-time dynamic 
approach we have to measure both the static and dynamic exponents together. It makes the 
work more difficult than measuring only the static exponents. However, it is clear if we are 
interested in both static and dynamic properties, the short-time dynamic approach is very 
efficient. The cluster algorithms are non-local and change already the dynamic universality 
class. On the other hand, the cluster algorithms can not straightforwardly apply to any 
systems, e.g., the systems with quenched randomness and lattice gauge theories. 



B. rUs^O 

Encouraged by the success for m^ = 0.0, we proceed to m^ ^ 0. The results in the 
last subsection show that the lattices sizes L = 64 and 128 are sufficient for our maximum 
updating times. The finite mo effect for mo = 0.01 is already invisible within the statistical 
errors. Therefore, in this subsection we will not systematically study the finite size and finite 
mo effects. However we have also performed some extra calculations to confirm the results 
presented. To locate the critical point Kc and estimate the exponent u, we need data for 
several K^s in the neighborhood of the critical point. The results in the last subsection also 
show that the difference of these K^s should be around half a percent of Kc or less. 

We perform simulations for rris = 0.2 and 0.4. For m^ = 0.2 we obtain data for K = 
0.2405, 0.2410, 0.2415 and 0.2420, which are shown in Fig. ^ with solid lines (from below) 
in log-log scale. For m, = 0.4 we obtain data for K = 0.3310, 0.3320, 0.3330 and 0.3340, 
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which are shown in Fig. ^ with sohd hnes (from below) in log- log scale. The lattice size is 
L = 128 and 2500 samples of the initial configurations are used for averaging. 

Following the procedure in the last subsection, searching for a curve with the best power 
law behavior in a time interval [tmic, 400], one can determine Kc and the critical exponent 6. 
Now we have data for four K^s and therefore have the choices to interpolate the staggered 
magnetization with three or four K^s. For Kc, both yield the same values within errors. 
With four K's, we obtain = 0.24153(12) for = 0.2 and = 0.33230(20) for = 0.4 
with tfnic = 20. The corresponding curves are shown in Figs. ^ and |^ with dashed lines. 
For rUs = 0.4, Kc is very stable for different choices of tmic- For = 0.2, Kc increases 
slightly when we take bigger t^ic- The reason is that our maximum updating times are not 
long enough. After carefully analyzing the data and looking at the curves in the figure, we 
convince ourself that Kc = 0.24153(12) is the most reasonable. 

In Fig. 1^, drM{t, T)\r=o is displayed with dashed lines for mg = 0.2 and 0.4 (from above). 
The results are obtained with four K^s. For = 0.4, the resulting exponent u is bad with 
three i^'s. The reason is clear. Kc is between K = 0.3320 and 0.3330. If we choose only 
three K^s, we have some systematic errors in interpolation. 

With Kc in hand, we perform simulations for mo = 0.0 and measure the second moment 
and the auto-correlation. A lattice size L = 64 is used and the number of samples of initial 
configuration for averaging is 7500. In Fig. |^ the second moment is displayed in log- log 
scale with solid lines for = 0.2 and = 0.4 (from above). For comparison, the second 
moment for = 0.0 is also plotted with a dashed line in the figure. Apparently, all three 
curves are parallel each other. The exponent y is independent of m^. In Fig. |^, the auto- 
correlation is displayed in log-log scale with solid lines for = 0.2 and = 0.4 (from 
below) and with a dashes line for rUs = 0.0. In all the relevant figures for non-zero m^, nice 
power law behavior is observed. 

In Table |, all the exponents are summarized. Within the statistical errors, the static 
exponent 2/5/z/ is independent of m^. The values also coincide well with those obtained 
from dynamics of model A and measured in equilibrium for the standard Ising model. The 
static exponent u for = 0.2 and 0.4 are 0.82(8) and 0.81(3) respectively. Even though 
the first value carries a relatively big error, they are clearly different from u = 0.64(2) for 
rUs = 0.0 and those measured from dynamics of model A and in equilibrium. However, 
u is independent of the value of the non-zero m^. These results are consistent with the 



theoretical calculations based on a theory coupling to a conserved current |^T[. With our 



accuracy, we can not detect any dependence of the dynamic exponent z on m^. This differs 
from the 0^ theory in Ref. [^. The exponent 6 is 0.129(5) and 0.148(4) for = 0.2 and 0.4 



respectively . These two values are obviously different from 6 = 0.108(5) for = 0.0, and 
show the nontrivial dynamic behavior for nonzero rris- However, with our data we should 
conclude that 6 depends on the value of the non-zero m^, even though there may be still 
some uncontrolled systematic errors. This point is different from the 0^ theory. 



IV. CONCLUDING REMARKS 

We have reported our Monte Carlo simulations of short-time critical dynamics for the 
three-dimensional anti-ferromagnetic Ising model with a globally conserved magnetization 
(not an order parameter). The power law scaling behavior has been observed for the stag- 
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gered magnetization (the order parameter), its second moment and the auto-correlation. 
All dynamic and static critical exponents are determined. For a conserved magnetization 
rris = 0.0, all the exponents are the same as those of model A. For a non-zero m^, the static 
exponents behave qualitatively the same as those of the 0^ theory coupling to a conserved 



current ||41[]. However, the dynamic exponents are somewhat different. The reason might be 
that the magnetization is not locally but only globally conserved. Therefore, it is interesting 
to study the dynamics, with a locally conserved order parameter or non-order parameter: 
in a flip, we exchange two spins in a local regime, the size of which must be much smaller 
than the lattice size. 

Acknowledgements: B.Z. thanks deeply Oerding for suggesting the topic and inspiring 
discussions. This work is supported in part by DFG; TR 300/3-1 and Schu 95/9-2. 
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FIGURES 
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FIG. 1. Time evolution of the staggered magnetization displayed in log-log scale. For solid lines, 
K is 0.22065, 0.22165 and 0.22265 (from below) for mo = 0.01, and K is 0.22115, 0.22165 and 
0.22215 for mo = 0.02. Dashed lines correspond to Kc = 0.22169(9) and 0.22163(5) for mo = 0.01 
and mo = 0.02. 
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FIG. 2. The second moment in log-log scale. For solid lines, K is 0.22065, 0.22165 and 0.22265 
(from below). The dashed line corresponds to Kc = 0.22164(22). 
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FIG. 4. The derivative drM{t,T)\T-=o in log-log scale. Solid lines are for nis = 0.0 and from 
below, mo = 0.01 and 0.02. Dashed lines are for = 0.2 and 0.4 (from above) with mo = 0.01. 
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FIG. 5. Time evolution of the staggered magnetization for = 0.2 in log-log scale . For 
solid lines, K is 0.2405, 0.2410, 0.2415 and 0.2420 (from below). The dashed line corresponds to 
= 0.24153(12). 
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FIG. 7. The second moment at Kc in log-log scale. For solid lines, rng = 0.2 and 0.4 (from 
above). The dashed line is for = 0.0. 
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TABLES 





Kc 9 X y Cd 


z 2/3/1/ u 


0.0 


0.22169(9) 0.108(5) 1.36(3) 0.965(11) 0.867(10) 


2.04(4) 1.03(4) 0.64(2) 


0.2 


0.24153(12) 0.129(5) 1.34(2) 0.954(4) 0.725(70) 


2.04(3) 1.05(3) 0.82(8) 


0.4 


0.33230(20) 0.148(4) 1.31(1) 0.945(8) 0.748(20) 


2.06(2) 1.05(3) 0.81(3) 


model A 


0.22170(4) 0.108(2) 1.36(2) 0.970(11) 


2.04(2) 1.034(4) 0.633(3) 


Ising 


0.221655(1) 


1.037(3) 0.630(1) 



TABLE L Critical exponents measured for different conserved magnetization m^'s. The values 
for dynamics of model A are taken from Ref. [46|. Under the item 'Ising' are the values for the 



standard Ising model in equilibrium |45| 
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